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The response of ultracold atomic Bose gases in time-dependent optical lattices is discussed based on direct 
simulations of the time-evolution of the many-body state in the framework of the Bose-Hubbard model. We 
focus on small- amplitude modulations of the lattice potential as implemented in several recent experiment and 
study different observables in the region of the first resonance in the Mott-insulator phase. In addition to the en- 
ergy transfer we investigate the quasimomentum structure of the system which is accessible via the matter- wave 
interference pattern after a prompt release. We identify characteristic correlations between the excitation fre- 
quency and the quasimomentum distribution and study their structure in the presence of a superlattice potential. 
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Ultracold atomic gases in optical lattices are a versatile 
laboratory for the study of fundamental quantum phenom- 
ena. The accurate control of the important physical param- 
eters over a wide range has been utilized for detailed exper- 
imental investigations of quantum phase transitions, e.g. the 
superfluid to Mott insulator phase transition (SF-MI) 
The primary experimental observable is the matter- wave inter- 
ference pattern of the atoms after release from the confining 
potentials and ballistic expansion. The interference pattern 
obtained after a sudden release of the atoms provides direct in- 
formation on the (quasi)momentum distribution of the system 
before the release. The interference pattern thus helps with 
the identification of different quantum phases, such as the SF 
and MI phases |0, . In recent experiments it was also used 
to study the response of the system in a two-photon Bragg- 
spectroscopy scheme based on a temporal modulation of the 
lattice potential. The broadening of the central interference 
peak after re-thermalization in a shallow lattice was used as 
a measure for the response . These techniques can also 
be employed in the context of more complicated irregular lat- 
tices as they can be produced, e.g., by using speckle patterns 
0,0, SHI] or two-color superlattice potentials HETIlll. In 
pioneering experiments on the response in two-color super- 
lattices the impact of a superlattice has been investigated 0]. 
In agreement with theoretical predictions a broadening of the 
response as function of modulation frequency was observed 
with increasing superlattice amplitude. 

Motivated by these experiments we study the response of 
a bosonic system to the modulation of the superlattice po- 
tential deep in the Mott-insulating regime. Going beyond 
the observables used in experiments, we study the quasimo- 
mentum structure of the system as it is revealed by the in- 
terference pattern after a prompt release of the atoms with- 
out any re-thermalization phase. To this end we perform an 
explicit time-evolution of the many-body state in the pres- 
ence of a time-dependent lattice based on the Bose-Hubbard 
Hamiltonian Similar studies have been done using the 
time-dependent density matrix renormalization group method 
(tDMRG) [13]. Our simulations reveal subtle correlations be- 
tween the modulation frequency within the resonance region 
and the quasimomentum distribution, which should be acces- 



sible to experiment. 

A system of N bosons in an one-dimensional superlattice 
potential with / sites at zero temperature is well described by 
the single-band Bose-Hubbard model (BHM) Qj]. The Bose- 
Hubbard Hamiltonian, formulated in second quantization in 
the basis of localized Wannier states of the lowest bands, reads 
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where the first sum runs over adjacent sites including a term 
connecting the first and last site of the lattice (cyclic boundary 
conditions). The Hamiltonian (Q]) with the bosonic creation 
(annihilation) operators k\ (a^) and the occupation number 
operator fi^ consists of three terms: The first term describes 
the tunneling between adjacent sites, the second term accounts 
for the on-site interaction of the atoms, and the third term rep- 
resents a site-dependent external potential. We introduce the 
superlattice potential via the latter term and describe its spa- 
tial structure by the reduced on-site energies q G [—1,0]. The 
physics of the BHM is governed by the competition between 
the relative strengths of these three terms, i.e. the tunneling 
strength J , the interaction strength /To, and the superlattice 
amplitude A . 

The basis of the BHM for N bosons and I sites is spanned 
by the occupation number states |{ni, 77,2, . • . , ni} a ) for all 
possible sets of occupation numbers rii with ^ rii = N. An 
arbitrary state can be expanded in this number- state basis lead- 
ing to |^) = J2a=i C <* K n i> n 2, • • • , ni} a ) with complex 
coefficients C a . The coefficients of the eigenstates \v) can 
be obtained by the numerical solution of the eigenvalue prob- 
lem of the Hamilton matrix. 

The basis dimension D grows factorially with / and N, 
thus limiting any calculation using the full basis to small sys- 
tems. However, in the strongly interacting regime (Uq ^> Jo) 
only a few basis states contribute to the low-lying eigen- 
states. This allows for a physically motivated truncation 
of the many-body basis. The relevant basis states 
are identified using the expectation value of the Hamiltonian 
({n 1 ,n 2 , • • • ,n/}| H | {ni,n 2 , • . • ,nj}) < E tmnc , and only 
states below the truncation energy E tmnc are included. By 
varying E tmnc one can explicitly assess and control the impact 
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FIG. 1: (color online) The lU-resonance in a system with / = N — 10 at Uo/Jo — 40 for a regular lattice (A = 0). Panel (a) illustrates the 
time- averaged energy transfer AE versus modulation frequency co of the lattice. Panel (b) shows the corresponding time-averaged interference 
pattern which also represents the quasimomentum distribution. Panel (c) shows the energy spectrum in the region of the first Hubbard-band 
with the vertical lines indicating sizable matrix elements (0| Hj | v) (darker shadings refer to larger values). The truncation energy of the basis 

is £ t runc/Jo = 120. 



of the truncation on observables. For regular lattices (Ao = 0) 
and filling factor N/ 1 = 1 this truncation allows for all rele- 
vant n-particle-n-hole states with n < E tmnc /Uo with respect 
to the reference state |l,l,...,l). 

We investigate the dynamics and response of the system 
induced by external time-dependent perturbations based on 
the explicit time evolution of the many-body state. Our cal- 
culations are motivated by recent experiments Q, [3] using 
a sinusoidal modulation of the lattice depth to perform two- 
photon Bragg spectroscopy. Unlike our calculations, these 
experiments include a rethermalization phase in the superfluid 
regime before the time-of-flight measurement in order to as- 
sess the energy transfer to the system. We assume a prompt 
release without rethermalization to directly access the quasi- 
momentum distribution after a certain modulation time. 

Formally, the time-dependent lattice potential is written as 
V(x, t) = [1 + F sm{ujt)] V(x), where V(x) is the static spa- 
tial lattice, uj the frequency, and F the relative amplitude of 
the modulation. All simulations are performed with a small 
amplitude F = 0.1 in accord with experiment. The time- 
dependence enters the Bose-Hubbard Hamiltonian H(t) via 
the time-dependent parameters J(t), U(t), and A(t), which 
are obtained within a Gaussian approximation for the local- 
ized Wannier functions H]. The parameters oscillate 
around their initial values Jo,Uq, and A at t = 0. 

We investigate the response of the system deep inside the 
Mott-insulating regime for fixed Uo/Jo = 40. The modu- 
lation frequency uj is varied in the range u/ Jo = 32 to 52, 
which covers the so-called 1U resonance at uo/ Jo ~ Uo/Jo — 
40. The ground state obtained for the initial Bose-Hubbard 
Hamiltonian Ho is used as initial state 0) and evolved up 
to t Jo = 10 for each frequency uo. The exact time-evolution 
is performed using a Crank-Nicholson scheme [01| Q with 
typically 3000 time steps. The relevant observables are evalu- 
ated at every 30th step. Following an initial phase with large 
changes, the observables saturate fl3ll and show only minor 
fluctuations Ill2ll . Since this residual time-dependence will not 
be resolved in experiment, we average the observables over an 
interval of evolution times within the saturated regime from 
t Jq = 6 to 10 (time-averaged quantities indicated by a bar). 



The simplest theoretical quantity to characterize the re- 
sponse of the system to the lattice modulation is the energy 
transfer AE(t) = (i/j,t\k \ ^,t)-(^,0|H | ^,0). As an ex- 
ample, Fig. QJa) shows the energy transfer for a system with 
/ = N = 10 in the Mott-insulator phase (U /Jo = 40) for 
a regular optical lattice (Ao = 0) in the vicinity of the 1U 
resonance at the modulation frequency uo = Uo- A detailed 
analysis of the energy transfer is given in Ref. 1121. 

Additional information can be obtained by a linear response 
analysis To this end, the Hamiltonian is lin- 

earized with respect to the modulation amplitude F, leading 
to Him(t) = Ho + FVo sin(o;£)[/iHo — kKj] with the am- 
plitude Vo of the lattice potential and the couplings fi and 
k. The first part of the linear term is irrelevant, because 
it only generates an energy shift. The second part couples 
the ground state to excited states via the tunneling operator 
Hj = —JJ2(ij) An a li near response picture, reso- 

nant transitions from the ground state |0) are expected when- 
ever the frequency uo coincides with the energy E v of an 
excited state \v) with a sizable matrix element (0| Hj| v). 
This interrelation is illustrated in Fig. QIc), where the energy 
eigenvalues for the first Hubbard band are shown. The verti- 
cal lines indicate excited states with sizable matrix elements 
(0| Hj I v). According to the linear response analysis, the en- 
ergies of those eigenstates should correspond to the resonance 
frequencies in the response. The comparison of the resonance 
energies emerging in the energy transfer in Fig. [Ha) with 
these excitation energies in[TJc) confirms this interpretation. 

In addition to the energy transfer we investigate the evolu- 
tion of the interference pattern. The matter- wave interference 
pattern after ballistic expansion is used to gain experimental 
information on the many-body state. In recent experiments 
0, [3] the central peak of the interference pattern was em- 
ployed as a measure for the response of the gas to a lattice 
modulation. The intensity distribution X(S) of the interfer- 
ence pattern as function of the relative phase 6 is given by 
2(8) = jEy^P^-i)^] (^| aja^l i/j) for an arbitrary 
state |^) 1 20]. For S = 2nq/I the intensity T{8) corresponds 
to the occupation numbers n q = (x/j \ c q c q | -0) =1(5 = 2irq/I) 
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FIG. 2: (color online) Interference pattern as function of uj in the 
region of the 1U resonance for a system with I — N = 10 at 
Uo/ Jo = 40 and Ao/ Jo = 2 for truncation energies E imnc / Jo = 120 
(a), E tmm / Jo = 80 (b), and E tmm / J = 40 (c). The color coding is 
the same as in Fig|4] 
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FIG. 3: Distribution of the a for / = 20 lattice sites of the incom- 
mensurate superlattice used in the calculations. 



of quasimomentum eigenstates with q = 5 I/2tt. Here, cj (c q ) 
are the bosonic creation (annihilation) operators with respect 
to the quasimomentum basis, which can be written as cj = 

TyE.expIi^aJ. 

Figure [TJb) illustrates the frequency-dependence of the in- 
terference pattern in the region of the 1U resonance. Each 
horizontal cut through the density plot represents the time- 
averaged interference pattern for a certain frequency uo. We 
assume that the lattice is switched off instantaneously after a 
certain evolution time in the modulated lattice — different from 
recent experiments which involve a re-thermalization period 
in the superfluid regime fl @]. The general interference struc- 
ture reveals a specific correlation between the frequency uo 
relative to the centroid of the 1U resonance and the quasimo- 
mentum distribution, i.e. the peaks of the interference pattern. 
Away from the resonance region, the intensity T(S) exhibits 
a broad background distribution characteristic for the Mott- 
insulating phase. For frequencies uo at the low-frequency end 
of the resonance a sharp interference peak emerges at S = 
indicating the resonant transition to the q = state. With 
increasing frequency uo this population moves to successively 
higher quasimomenta \q\, i.e. the interference peak splits and 
shifts towards larger \S\. The fine- structure of the resonance 
is thus mapped onto the interference pattern in an experimen- 
tally accessible way. 

So far, the simulations were restricted to small systems with 
large bases. To treat larger lattices we reduce the trunca- 
tion energy i^nmc controlling the basis size. Figure [2] illus- 
trates the insensitivity of the interference pattern on changes 
of i^tnmc- There is practically no difference when reducing 
the truncation energy from E tmnc /Jo = 120 to 80. Even 
for E tmnc / Jq = 40 all relevant features are reproduced, al- 



though the intensity of the interference peaks shows slight de- 
viations. All qualitative conclusions regarding the correlations 
between frequency and quasimomentum distribution remain 
unaffected. 

Using the lowest truncation energy we investigate the re- 
sponse of a system with I = N = 20 at Uq/Jq = 40. The major 
focus is on the change of the response and the interference 
pattern if a two-color superlattice potential of increasing am- 
plitude Ao is added. The distribution of the relative strengths 
€i for the incommensurate superlattice used in the following 
are shown in Fig. [3] The the ratio of the wavelengths of the 
two standing waves isAi/A2~0.81, similar to a recent exper- 
iment 10]. 

Figure [4] depicts the evolution of the interference struc- 
ture as function of the superlattice amplitude Ao for the in- 
commensurate case. The right-hand panels show the energy 
spectrum with vertical lines marking sizeable matrix elements 
(0| Hj| v). The left-hand panels depict the energy transfer 
as function of frequency. The result for Ao/ Jo = shown 
in Fig. |4ja) confirms our previous discussion of the smaller 
regular lattice in Fig. [T] The linear response analysis pro- 
vides a good estimate for the resonance energies via the ener- 
gies of those excited states that exhibit strong transition ma- 
trix elements to the ground state. Furthermore, the correlation 
between frequency uo and the quasimomentum distribution is 
even more pronounced. At the low-frequency end of the res- 
onance quasimomentum states around q « are populated, 
whereas for larger modulation frequencies successively higher 
quasi-momenta are occupied. 

The results for a small superlattice amplitude Ao/ Jo = 1 in 
Fig. Hfb) show minor changes in the quasimomentum struc- 
ture and the energy transfer as compared to (a). Nevertheless, 
the number of possible excitations from the ground state in- 
creases. A further increase of the superlattice amplitude to 
Ao/ Jo = 2 leads to a weak suppression of the interference 
structure as shown in Fig. HJc). In comparison to the energy 
spectrum in (b) there are many more large matrix elements 
which are not localized at distinct energies but spread over the 
whole range. The occurrence of the small gaps at both ends 
of the energy band is also visible in the density plots as small 
shifts in the interference structure along the energy (vertical) 
axis. This also indicates the broadening of the resonance due 
to the superlattice lB[[l21. 

Further increase of the superlattice amplitude to Ao/ Jo = 4 
leads to the disappearance of the interference peaks as de- 
picted in Fig. Hfd). The number of strong matrix elements 
which couple ground and excited states is further increased. 
Although the interference pattern does not show sharp peaks 
anymore, the energy transfer in these cases still exhibits a 
strong resonance behavior Jl2ll . Only the fine-structure of the 
resonance in the energy transfer is affected by the superlat- 
tice. However, the distinct correlations between excitation fre- 
quency and quasimomentum distribution vanish far below the 
transition from the homogeneous Mott insulator to the Bose- 
glass phase at Ao ~ Uq. These general features are robust 
against changes of the superlattice structure. We have per- 
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FIG. 4: (color online) Energy transfer (left-hand panels), interference 
pattern (middle panels), and excitation spectrum (right-hand panels) 
of a system with I=N=20 and Uo/ Jo =40 for several superlattice 
amplitudes Ao as indicated in the plots. The density plots illustrate 
the correlations between quasimomentum distribution and excitation 
frequency in the region of the lU-resonance. The vertical lines in the 
energy spectra (right-hand panels) point out strong matrix elements 
<0|Hj|i/>. 



formed explicit simulations using a commensurate superlat- 
tice showing the same behavior. 

The rapid change of the interference pattern can be ex- 
plained in the linear response picture: In the absence of a 
superlattice, resonant transitions connect the ground state to 



a few excited states only — those characterized by large tran- 
sition matrix elements (o| Hj| v) [Fig. 4(a) and (b)]. The 
many-body state \ip,t) during the time evolution is domi- 
nated by these few states and exhibits well-defined interfer- 
ence peaks. With increasing superlattice amplitude Ao more 
and more sizable transition matrix elements emerge and the 
time-evolved state is a superposition of a large number of 
eigenstates [Fig. 4(c) and (d)]. The fragmentation of the 
state causes a fragmentation of the interference pattern and 
in effect a suppression of the distinct peaks. This mechanism 
can be confirmed within a simple toy-model by comparing the 
interference pattern of an excited eigenstate with the one re- 
sulting from a coherent superposition of a few neighboring 
eigenstates. Note that this phenomenon is quite different from 
the Mott-insulator to Bose-glass transition, which appears at 
much larger superlattice amplitudes. 

In summary, we have discussed the response of Bose gases 
in modulated lattice and superlattice potentials with an em- 
phasize on the quasimomentum distribution and the interfer- 
ence pattern after prompt release from the lattice. For mod- 
ulation frequencies in the region of the 1U resonance distinct 
peaks appear in the interference pattern. Their position is cor- 
related with the modulation frequency relative to the centroid 
of the resonance: with increasing frequency they shift to larger 
quasimomenta. In the presence of a superlattice this character- 
istic correlation vanishes already for small superlattice ampli- 
tudes, much faster than the resonance observed in the energy 
transfer. 
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